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The impact of mergers on relaxed X-ray clusters 
III. Effects on compact cool cores 
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ABSTRACT 

We use the simulations presented in lPoole et al.l (|2006f ) to examine the effects of merg- 
ers on the properties of cool cores in X-ray clusters. Motivated by recent Chandra and 
XMM observations, we propose a scheme for classifying the morphology of clusters 
based on their surface brightness and entropy profiles. Three dominant morphologies 
emerge: two hosting compact cores and central temperatures which are cool (CCC 
systems) or warm (CWC systems) and one hosting extended cores which are warm 
(EWC systems). In the cases we have studied, CCC states are disrupted only after 
direct collisions with merging cluster cores. This can happen in head-on collisions or 
during second pericentric passage in off-axis mergers. By the time they are relaxed, 
our remnant cores have generally been heated to warm core (CWC or EWC) states but 
subsequently recover CCC states by the end of the simulation. The only case resulting 
in a long-lived EWC state is a slightly off-axis 3:1 merger for which the majority of 
shock heating occurs during the accretion of a low-entropy stream formed from the 
disruption of the secondary's cool core. Since td yn <C t coo i for all our relaxing merger 
remnant cores, compression prevents their core temperatures from falling until after 
they relax to the compact states allowed by their remnant central entropies. This nat- 
urally explains the population of relaxed CWC systems observed in recent Chandra 
and XMM observations with no need to invoke AGN feedback. T he morphological 
segreg ation in the L x — T x and f3 — r c scaling relations noted by iMcCarthv et alA 
(2004J) are qualitatively reflected in the results of our mergers as well. However, none 
of the cases we have studied produce systems with sufficiently high central entropies 
to account for the most under- luminous EWC systems observed. Lastly, mergers do 
not efficiently mix the ICM. As a result, merging systems which initially host central 
metallicity gradients do not yield merger remnants with flat metallicity profiles. Taken 
together, these results suggest that once formed, compact core systems are remarkably 
stable against disruption from mergers. It remains to be demonstrated exactly how 
the sizable observed population of extended core systems was formed. 

Key words: cosmology: theory - galaxies: clusters: general - intergalactic medium 
- X-rays: general 



1 INTRODUCTION 

The latest generation of X-ray telescopes have revolution- 
ized our understanding of the intracluster medium (ICM) of 
galaxy clusters but yet one of the oldest observations of the 
field remains poorly understood. As significant sample sizes 
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of X-ray clusters initially became available, it was quickly 
realized that, broadly speaking, two populations could be 
identified fro m the shape of a system's X-ray surface bright - 
ness profile l| Jones fc Formanl ll98J : lOta fe Mitsudall2004f ): 
one with highly peaked central surface brightnesses (com- 
pact cores; r c ~ 50kpc) and one with flat central surface 
brightnesses (extended cores; r c ~ 150 kpc). It was quickly 
understood that high and centrally peaked gas densities are 
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primarily responsible for the central surface brightnesses of 
compact cores while extended cores possess lower and more 
constant central gas densities. 

As a consequence of the p 2 dependence of the ICM's 
emissivity, radiative cooling should be significantly more rel- 
evant to the structure of compact cores. This assertion is 
supported by observations which have indicated that com- 
pact co res tend to host centr al positive temperature gra- 
dients (|De Grandi fc Molendil 120021 ). low central entropies 
^McCarthy et al. 112004 ') and radiative cooling rates which 
should release their thermal energy on timescales (the "cool- 
ing time "; t coo i herea fter) much shorter than the mean age of 
clusters (|Allenlll998l ). Within the context of specific physical 
models, these systems were once thought to host bulk radial 
flows of cold gas arising from the loss of pressure support 
their short cooling timescales induce. As a result, they were 
classically referred to as "cooling- flow" systems. Extended 
cores on the other hand tend to have flat central temper- 
ature profiles, elevated central entropies and cooling times 
in excess of a Hubble time. Classically, these systems were 
referred to as "non-cooling flow" systems. 

Great observational progress has been made since the 
launch of Chandra and XMM and cluster cores are now 
understood to be much more complicated than previously 
believed. For example, it has been shown that clusters with 
compact cores can in fact have a range of central tem- 
peratures - some significantly cooler than their surround- 
ings and some rough l y isothermal (|Donahue et al] 1200a : 
ISanderson et al] 120061 ; IPratt et all 120071 ). Hence, compact 
cores are not necessarily cool. It is presently unclear what 
relative roles active galactic nuclei (AGN), mergers and cool- 
ing play in setting this spectrum of states. 

As a result, the terminology of the field is in flux. Given 
this situation, we propose here the adoption of new nomen- 
clature for describing the morphology of X-ray cluster cores. 
Anticipating the results to be presented in this paper, we 
will refer to extended core clusters as extended warm core 
(EWC) systems and compact core clusters as either com- 
pact cool core (CCC) or compact warm core (CWC) sys- 
tems. We shall illustrate in this paper why this facilitates a 
more insightful connection to the physical processes shaping 
the structure of cluster cores. 

Despite recent progress in understanding the structure 
of cluster cores, the origin and precise nature of their ob- 
served distribution of surface brightness core radii still lacks 
a clear and definitive explanation. Several hypotheses have 
been proposed to account for these populations, with three 
being particularly persuasive. 

First, it has been known for a decade that radiative 
cooling and shock heating during accretion are not suffi- 
cient to_accoun^or_theobserved structure of X-ray clus- 
ters. iMcCarthv et al. I (|2004 M04 henceforth) have devel- 
oped a model in which entropy injection (possibly due to 
galactic outflows or early AGN activity) during the proto- 
cluster phase establishes a spectrum of minimum entropies 
for clusters leading to variations in cluster core cooling ef- 
ficiencies. Reasonable amounts of such additional entropy 
(~ 200 keVcm 2 ) and radiative cooling can account for the 
normalization and scatter in the observed L x — T x and 
Mt — T x relations while obeying observed cold gas con- 
straints. M04 also illustrate how clusters receiving little 
entropy injection quickly evolve to states resembling com- 



pact cool core systems as a natural consequence of radia- 
tive cooling. Systems receiving higher levels remain in long 
lived extended core states. Although not stated explicitly 
in their discussion, systems evolving to compact cool core 
states in this model would then be stabilized through sev- 
eral possible feedback processes involving late AGN activity 
JVoit fc Donahuell200a ; iNusser et aill2006l : IMcCarthv et all 
12007(1 

Second, the entropy injection required to account for 
the structure of extended cores may have occurred more 
recently through interactions of AGN with the ICM. Such 
interactions have been observed in the cores of several 
systems where jets from AGN are seen to be actively inflat- 
ing large high-entropy cavities (or "bubbles"). It has been 
shown that the power supplied by such jets correlates with 
the local Bondi accretion rate and is energetically capable 
of offsetti ng the present-day effects of cooling (lAllen et al 
| 2004 iBenson et all 120031; ICroton et all 1200a : iBower et°af 



12009 : JNemmen et alJuOOTfl Typically, such interactions are 
localized to the central ~ 50 kpc of a system but in a few 
instances, much larger outbursts reaching to ~ 200 kpc have 
been o bserved ( McNamara ct al. 2005). However, as pointed 
out by IMcCarthv et all l|2007T ). present-day AGN activity is 
energetically capable of merely maintaining the present day 
configuration of actively cooling cluster cores. The energy 
required to transform such systems into the observed range 
of non-cooling cores is much greater than the largest AGN 
outbursts yet observed. This fact strongly suggests that 
an early epoch of " pre-heating" due to a yet unspecified 
mechanism (possibly AGN) is largely responsible for setting 
the observed range of present-day cluster core states. 
Unfortunately, the short timescales, complicated magneto- 
hydrodynamic physics and large range of relevant physical 
scales involved in the interaction of jets with the ICM makes 
exploring su ch issues from theoret i cal directions extremely 
challenging (IChurazov et all 2001 : Rus zkowski et al. 2004; 
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120051; IVernaleo fc Reynolds! I2006T ; iHeinz et all 120061 ; 
battaneo fe Tevssierl 120071 ). 

lastly, mergers between clusters offer a compelling 
means by which to transform compact core systems into 
extended core systems. It has been a tacit belief for some 
time that extended core systems are a product of merger 
interactions but the issue has received remarkably little 
detailed consideration. This hypothesis is primarily based 
on early anecdotal observations that few compact cores 
exhibited disturbed morphologies while the limited num- 
ber of studied extended core systems tended to be dis - 
turbed (|Edge. Stewart, fe Fabianlll992l ; lBuote fe Tsailll996l ). 
However, with improved instrumentation and larger sample 
sizes, there are now several examples of d i sturb ed compact 
core systems {e.g. A2204, ISanders et all 120051 ) . including 
ones presently undergoi ng significant mergers {e.g. A2142, 
iMarkevitch et all 120001 ). as well as extended core sys- 
tems lacking any evidence of having been disturbed at all 
{e.g. Abell 2034 and Abell 2631; K. Cavagnolo, private com- 
munication) . 

Despite such concerns, the merger-origin hypothesis for 
extended core systems is appealing on theoretical grounds. 
The presently favored hierarchical clustering paradigm for 
structure formation predicts that large objects are formed 
through regular accretion of smaller subclusters with each 
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observed system representing the product of a unique his- 
tory of many such events. This process is seen to be active 
to the present epoch with recent Chandra and XMM ob- 
servations routinely detecting evidence of significant recent 
or ongoing mergers. With as much as 10 64 ergs of initial 
kinetic energy, mergers are energetically capable of signifi- 
cantly modifying the structure of a system's ICM through 
the shocks they produce. Furthermore, the wide variety of 
individual accretion histories may lead naturally to a range 
of remnant states like that described by M04. 

Recent studies that have investigated cluster mergers 
in a context capable of testing the merit of this scenario 
have yielded inco n clusiv e results. In the influential study of 
iRicker fc Sarazinl i|200ll . RS01 hereafter) for instance, some 
evidence for the flattening of central density profiles due 
to shock heating from mergers is presented. However, since 
the focus of their study was on merger induced luminos- 
ity and temperature "boosts", they did not include the ef- 
fects of cooling (which are known to be important to the 
structure of compact core systems) and they present rela- 
tively little detailed analys is of their remnant cores ' struc- 
tures. iGomez et aT\ (|2002l ) and iRitchie fc Thomas I (|2002h 
have constructed suites of idealized merging systems but 
their conclusions regarding the st ability of co mpact core sys- 
tems to mergers seem to conflict. IGomez et al. I (| 20021 ) find 
that systems with initially sho rt cooling times wil l likely 
reestablish their structure while IRitchie fc Thomas I (J2002J) 
find that compact core systems are completely disrupted in 
equal mass and near-axis 8:1 mergers. Adding to the con- 
troversy, studies of clusters produced in cosmological hydro- 
dynamic simulations routinely fail to yield enough scatter 
in th eir X-ray scaling re l ations to accoun t for observations 
(e.g. iRowlev et ai1l2004l : iKav et alj|2007l) . 

To address this and other gaps in our understanding 
of cluster mergers, we have initiated a numerical study of 
their impact on the observable properties of relaxed cool 
core X-ray cluster s using idealized tw o-body simulations. In 
a previous paper IjPoole et alj|200a , Poole06 hereafter) we 
have presented the details of our numerical approach and 
the dynamical evolution of a suite of 9 simulations (three 
mass ratios with three impact parameters each). We then 
used these simulations for a second study examining the 
effects of mergers on global X-ray and SZ observables and 
their scaling relations (jPoole et al.ll2007l . Poole07 hereafter) . 
The same suite of simulations form the basis of the study 
presented in this paper. We briefly review several aspects of 
our simulations particularly relevant to this paper in Section 
[5] and refer the reader to Poole06 for more details. 

In this study we establish criterion for quantifying the 
evolving morphology of merging systems (Section 13. 2[) and 
present their evolution through three distinct morphological 
states (Section 13. 4p . We then relate the occurrence of these 
states to changes in the structure of the system's core (Sec- 
tion 13. 5[) . After considering the morphological segregation 
mergers induce in observed scaling relations (Section 13, 6p 
we then consider the mixing effects of mergers on metallicity 
gradients (Section [4}. We then discuss several consequences 
of our findings in Section [S] and finally, summarize our study 
in Section [6] 

In all cases our assumed cosmology will be 
(n A /,O\) = (0-3,0.7) with H =75kms- 1 Mpc~ 1 . 



2 SIMULATIONS 

We run our simu l ations with GASOLINE 
(|Wadslev. Stadel. fc Quinnl [2004), a versatile parallel 
SPH tree-code with multi-stepping. We include the effects 
of radiative cooling in our simulations and feedback from 
star formation but the effects of jets from active galactic 
nuclei (AGN) are omitted. 

The basis of our study is a set of 9 idealized 2-body 
cluster merger simulations initialized to accurately match 
the properties of observed compact cool core clusters. Three 
mass ratios are examined (1:1, 3:1 and 10:1) with the most 
massive system set to have a mass of 1O 15 M0 in every case. 
For each mass ratio, three impact parameters are examined 
(one head-on and two off-axis cases). We parametrize impact 
parameters by the smaller system's transverse velocity when 
its core crosses the virial radius of the massive system, in 
units of the circular velocity of the massive system at that 
radius (vt/V c ). 

In Poole06 we present a detailed account of how our 
isolated clusters and their orbits are initialized as well as the 
numerical methods and code parameters of our simulations. 
In this section we summarize key aspects of our approach 
pertinent to our present analysis and direct readers looking 
for additional details to Poole06. 



2.1 Initial conditions 

To initialize the structure of our s ystems, we h a ve fol - 
lowe d the analytic p r escrip tion of iBabul et al. I (| 20031 ) 
and iMcCarthv et al. I (J2004) to produce cool core clus- 
ters which conform with recent theoretical and observa- 
tional insights into cluster structure. The dark matter 
density profiles of our systems follow an NFW- like form 
|Navarro. Frenk. fc Wh"it3ll996l : lMoore et aUH998h with the 
central asymptotic logarithmic slope chosen to be = 1.4 
and the scale radius (r s ) selected to yield a concentration 
c = R20o/r s — 2.6 (we will use i?A throughout to indicate 
the radius within which the mean density of the system is A 
times the critical density, p c = 3Hq/8tvG; R200 — 1785 kpc, 
^500 = 1166 kpc and R coo i — 180 kpc initially for our most 
massive systems). The initial gas density and temperature 
profiles of the clusters are set by requiring that (1) the gas 
be in hydrostatic equilibrium within the halo, (2) the ratio 
of gas mass to dark matter mass within the virial radius 
be f2f,/(f2 m — fib), and (3) the initial gas entropy scale as 
S(r) oc r 1 * over the bulk of the cluster body, as found in 
cosmologi cal hydrodynamic simulations and observations of 
the I CM (|Lewis et al.ll2000l : IVoit et al.ll2005l : TPonahue et ail 
120061 ). We normalize the entropy profiles such that the tem- 
perature of the ICM at R V i r is half the virial temperature. 
We construct orbits for our systems which produce spec- 
ified radial and tangential velocities for the secondary sys- 
tem (v r and vt respectively) when its centre of mass reaches 
the virial radius of the primary (R V i r ). For each of the three 
mass ratios we study, we examine three orbits selected to 



6 We use the standard proxy for entropy given by S = kT/nJ 
with n e and T representing the electron density and temperature 
of the gas. 
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produce a typical value of v r (Rvir) and to cover a signif- 
icant range of the transverse velocity Vt(Rvir) giving rise 
to mergers found in cosmological d ark matter simulations 
|Tormenlll997l ; IVitvitska et alil2002l) . 

Throughout our analysis we shall distinguish the two 
merging systems by calling the more massive system the 
"primary" system (we arbitrarily choose one to be the pri- 
mary in the 1:1 cases, with no consequences for our results 
given their high degree of symmetry) and the less massive 
system as the "secondary" . 



2.2 Dynamical evolution 

In Poole06 we illustrated the generic dynamical progression 
which each of our simulations proceed through, identifying 
five distinct stages: a pre-interaction phase, first core-core in- 
teraction, apocentric passage, secondary core accretion, and 
relaxation. These stages feature prominently in the discus- 
sion which follows so we will briefly review their progression 
here, for those who have not read Poole06. 

As the cores accelerate towards each other during pre- 
interaction, a pair of shock fronts materialize and are driven 
towards each core, heating and compressing them briefly. 
The core of the secondary system crosses within -R200 of the 
primary system at t and at t c i osest they reach closest ap- 
proach (i.e. first pericentric passage). In every case (includ- 
ing the head-on collisions) , some part of the secondary's core 
survives its first encounter with the primary, tidally stripped 
into large cool streamers in off-axis cases and strings of sev- 
eral clumps in the on-axis cases. At t apo the disturbed cores 
reach maximum separation. After subsequently reaching sec- 
ond pericentric passage at t accre t e , no observable trace of the 
secondary core remains. This period marks the beginning of 
a prolonged period of accretion when the plume of mate- 
rial generated from the disruption of the secondary's core 
accretes as a high velocity stream. This instigates a second 
episode of core heating followed by a period of relaxation. 
At approximately t re i a x, the system exhibits no obvious sub- 
structure in simulated 50ks Chandra images with the system 
placed at z — 0.1 (see Poole06 for sample images and details 
on how this is done) . 



2.3 Consequences of our coarsely sampled initial 
conditions 

In this paper, we seek to understand the role which merg- 
ers are playing in shaping the observed diversity of clus- 
ter core morphologies. It is important to note however that 
our simulations are certainly not uniformly represented in 
the statistics of the population. In Poole07 we used merger 
rate statistics from cosmological N-body simulations to ar- 
gue that our 1:1 mergers are expected to be rare, our 3:1 
mergers uncommon but of statistical significance and our 
10:1 mergers extremely common (see Poole07 for a more 
quantitative account). For these reasons, when considering 
the effects of mergers in the context of statistical properties 
of cluster populations, we will be primarily concerned with 
the effects of 3:1 and 10:1 events. 

We will find in later sections that the structure of clus- 
ter cores following mergers is the result of a delicate balance 
between shock heating and radiative cooling. In particular, 



we shall illustrate in Section [3.5.1 1 the important role stream 
accretion (see Poole06) plays in setting the entropy of our 
merger remnant's cores. The efficiency of this process de- 
pends sensitively on many factors. These include the mass 
and density of material present in the stream and the veloc- 
ity it carries to the remnant core. Both of these depend on 
the initial mass and orbit of the secondary in a complicated 
way, through (for example) their influence on the secondary 
core's disruption and the amount of dynamical friction ex- 
perienced by the system. 

The impact parameters and mass ratios of our simula- 
tions were selected to represent the most interesting range of 
situations exhibited by dark matter substructure in cosmo- 
logical simulations. However, these ranges are sampled very 
sparsely by our simulations. As a result, we do not claim 
to have performed a systematic analysis of many issues dis- 
cussed in this paper. We shall merely seek to reveal the gen- 
eral processes most relevant to establishing the structure of 
cluster cores following mergers. Many details will have to be 
left open to future study. 

2.4 A cautionary note 

In SPH cluster simulations, cooling converts hot gas in the 
centre of the system into dense cold gas which is then (if 
star formation is included) transformed into a collisionless 
stellar medium. As a result of the finite timescales involved 
in transforming cold gas into stars, a multiphase interface 
forms between dense cold gas and hot rarefied gas at the cen- 
ter. It is well known that multiphase in terfaces of this sort 
are poorly treated in SPH simulations (JRitchie fc Thomas I 
l200ll : lMarri fc Whitell2003r ). Since the dynamical time in the 
center is much shorter than the cooling time, the result is a 
quasi-adiabatic flow which drives compressional heating as 
hot material collects several kernel radii outside of the mul- 
tiphase interface. This leads to artificially suppressed densi- 
ties and enhanced temperatures and entropies in the cores. 
The pooling of baryons to the center also leads to deeper 
central potentials which contribute to higher core tempera- 
tures. During steady undisturbed cooling in our simulations, 
these effects appear initially at r ~ 3kpc but slowly prop- 
agate outwards with time. 

Such numerical concerns are significant only while a 
core is dense and its cooling timescales very short. In each 
of our merger simulations, the cores of both systems become 
sufficiently disturbed by the interaction to restrict periods of 
active quiescent cooling to durations shorter than 4 Gyrs. To 
determine the physical extent of the regions influenced by 
this artifact, we have run our systems in isolation at both our 
nominal resolution and with 4x its mass resolution. After 
4 Gyrs, we find that the density and temperature profiles of 
these two simulations agree to within 10% outside of 30kpc. 
Thus, the influence of such numerical effects should remain 
confined to radii < 30 kpc, marginalizing these concerns for 
our study. 

To ensure that concerns which might arise as a result 
of these effects are minimized, we will exclude the central 
30 kpc when studying projected quantities (such as surface 
brightnesses and central temperatures) and measure "cen- 
tral" quantities (such as entropy) outside this radius, at 
40 kpc. At a redshift of z — 0.1, 30 kpc corresponds to 17"; 
comparable to the typical central bin sizes of Chandra tem- 
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perature profiles or the ROSAT PSF (~ 15"). Since these 
are the instruments used to obtain the data against which 
we compare our results, removal of this region should not 
introduce significant biases during our comparisons to ob- 
servations. 

Since the effect of rising central temperatures is to re- 
duce central temperature gradients, timescales of the recov- 
ery of warm core to cool core morphologies are likely exag- 
gerated in our simulations, leading to systematic biases in 
our analysis. It is critical to note however that the direction 
of this bias is always the same: towards overestimating the 
longevity of warm core morphologies. Hence, the broad un- 
derstanding and conclusions which emerge from this study 
are preserved. 



3 ANALYSIS 

In this section we present our methods of classifying the 
evolving morphologies of our simulated mergers. We shall 
find that as our mergers progress, they exist primarily in 
one of three states: compact cool core (CCC), compact warm 
core (CWC) and extended warm core (EWC) morphologies. 
In the sections which then follow, we shall discuss the tem- 
poral evolution of our mergers through these three states 
and relate their occurrence to the changing structure of our 
merger remnants' cores. 



3.2 Classification of Morphological States 

We seek to determine how our evolving merger systems com- 
pare with recent Chandra and XMM observations of the 
ICM: at which points they would manifest as compact or ex- 
tended core systems, when they would appear to have been 
significantly affected by radiative cooling based on their en- 
tropy and projected temperature profiles, and when they 
would exhibit evidence of merger activity under typical ob- 
servational circumstances. To minimize biases in our com- 
parison to observations, we have developed a classification 
scheme based on a set of easily observable properties of the 
system which can be readily related to its physical state. 

It has been noted that projection effects can have 
a significant influence on the surface brightness struc- 
ture and projected tem peratures of a merging system 
IjRitchie fc Thomas 1120021 ). For this reason, we have consid- 
ered three orthogonal projections in our study: two in the 
plane of the orbit along the line initially connecting the sys- 
tems (x) and in the direction of initial tangential velocity 
(y) and one normal to the plane of the orbit (z). However, 
given the rarity of observing systems in a configuration cor- 
responding to our x-projections, we shall focus only on the 
y and z-projections in this paper. The only significant effect 
of this omission is to exclude those rare and peculiar states 
produced when the systems are separate and still relaxed 
but appear as one system in projection. 



3.1 Observational Datasets 

To compare our classified systems to observations, we 
will employ the 20 clu ster statistical sample presented by 
ISanderson et alj (120061 ) . We have also augmented the results 
of this study by fitting 1 dimensional j3- models (c.f. Equa- 
tion 1 below) to each cluster. To do so, we took the APEC 
spectral normalization per unit area as a measure of pro- 
jected emissivity, using th e more finely binned, p rojected 
annular spectral profiles of ISanderson et al.l (|2006f ) to con- 
struct azimuth ally-averaged radial su rface brightness pro- 
files. Eqn. 1 of ISanderson et al.l (|2006f ) was fitted using the 
X 2 statistic to each profile across the full radial range to 
determine the core radius; no attempt was made to excise 
or model any excess emission associated with a cool core. 
The use of projected emissivity instead of actual X-ray sur- 
face brightness is advantageous, since it allows for the ef- 
fects of galactic absorption, which is a free-fitted parame- 
ter in each annular spectrum for most of the clusters ana- 
lyzed (but fixed at the level expected from HI surveys in the 
few clusters where this pa rameter was not well constrained; 
see ISanderson et al.l 120061 . for details). For example, Abell 
478, suffers from strong galactic absorption with significant 
spatial variation (Sande rson et alj |2005) that could bias fits 
based on the surface brightness. 

In later sections we also employ th e catalog u e of R OSAT 
and ASCA observations presented by iHorner I (|200ll ). This 
catalogue is best suited to comparisons of the global X-ray 
properties of our simulations to observations because it is 
large, complete and presents results in which the cores of 
systems are not excised (a common procedure in other cat- 
alogues performed to reduce scatter in scaling relations for 
cosmological studies). 



3.2.1 surface brightness 

We first classify our systems according to the shape of their 
projected surface brightness profiles. To do so, we utilize the 
simple 3-parameter /3-model commonly employed in analysis 
of observations. This model characterizes the azimuthally 
averaged surface brightness profile of a system at a projected 
radius 7 r p by the equation 



2-fx\Tp) — 



E 



(l + (r P /r c rr- 1)/2 



(1) 



where r c is the profile's core radius, /3 characterizes its 
asymptotic logarithmic slope at large radius and E D its nor- 
malization. We exclude projected radii less than 30 kpc from 
this analysis and have verified that varying this radius from 
20 to 50 kpc has a minimal effect and does not influence the 
conclusions of our study. 

Although more sophisticated parameterizations of the 
X-ray surface brightness profiles of clusters have been de - 
veloped (e.g. IVikhlinin et al.l 120061 : [Mahdavi et HI 120071 ). 
we utilize this one to take advantage of the large statistical 
samples of /3-model fits presently available. 

To separate "compact" cores from "extended" cores we 
shall use a cut of 100 kpc. This corresponds approximately 
to the unp opulated region of the observed r c distribution re- 
ported by lOta fc Mitsudal l|2004h separating observed com- 
pact and extended core systems. 



7 Throughout this paper, we adopt the peak of the projected 
bolometric X-ray surface brightness as the center of the system. 



6 Gregory B. Poole et al. 



3.2.2 central entropy and temperature gradient 

To judge the relevance of radiative cooling for our evolv- 
ing cluster cores we will use the central entropy of our sys- 
tems. We do this by computing radially binned entropy pro- 
files, denoted as S(r), and interpolate from them to obtain 
"central" values at 40kpc (denoted as S40). We use entropy 
for this purpose rather than the system's projected tem- 
perature because it is numerically robust in our simulations 
and unambiguously tied to the relative importance of heat- 
ing and cooling for both our observed and simulated sys- 
tems. Furth ermore, as pointed out p reviously by several au 



thors 



£R 



rmore, as pointed out p reviously by several au- 
Voit et aliliool. 120031 : iKaiser fc Binnevl 120031 ; 



iMcCarthv et al.l 120071 ). variations in the core properties of 
X-ray clusters are most insightfully understood in terms of 
this quantity. 

To separate "cool" core systems from "warm" core sys- 
tems, we shall use a cut in S40 of 50keVcm 2 and refer to 
systems with S40 < 50keVcm 2 as "cool core" systems and 
systems with S40 > 50keVcm 2 as "warm core" systems. 
This choice is motivated by an exam ination of the entropy 
profiles of the ISanderson et al.l (|2006r ) catalogue (these will 
be presented in a future publication) ; central entropies lower 
than this value tend to occur in cool core systems while 
higher central entropies tend to occur in warm cores. Fur- 
thermore, recent examinations of BCGs in galaxy clusters 
find a break at this entropy level such that those in clusters 

with £40 < 50 keVcm 2 show evidence of star formation while 

~ 1 1 

those in systems with wa rmer cores do not (Bildfell et al.l 

120081 : iRaffertv et alJlJOOSh . For reference, the cooling time 

of material with S = 50 keVcm 2 is ~ 2 Gyrs for our systems. 

However, temperature is an intuitively appealing quan- 
tity and can be more directly compared to observations. 
For this reason we have also measured the shape of our 
systems' evolving projected central temperature profiles. 
Chandra and XMM observations have shown that pro- 
jected cluster temperature profiles have a universal shape 
outside of the projected radius O.I-R200; a broad peak be- 
tween 0. 1-R2 00 and 0.27^200 and a continuous decline outside 
of O.2.B200 ( Markevitch et al.lll99Sl ; ISanderson et al.1 l200rj ; 
IPratt et alj 120071 ). For projected radii within O.I-R200, a 
great deal of system-to-system variation in temperature pro- 
files is seen. This radius roughly corresponds to the cooling 
radius of compact core clusters and it is within this region 
that cool core system s exhibit relatively low temperatures 
jPiffaretti et al.ll2005h . 

Motivated by these observations, we have chosen to 
characterize the apparent temperature structure of the cores 
of our merging systems by measuring the strength of their 
central temperature gradients with the quantity 



A T = 



J- o u 



J- it 



-L 011 



(2) 



Here, Ti nner is the spectrally fit temperature (see Poole07 
for details regarding our method of simulating observed clus- 
ter temperatures) of material in projected radii from 30 kpc 
to O.I-R200 (0 to O.I-R200 f°r the observations; this repre- 
sents a negligible difference however) and T ou t er the spec- 
trally fit temperature for material in the range of projected 
radii O.I-R200 to 0.27?2oo- This quantity is thus a measure 
of the temperature gradient within the core region, where 



system-to-system variations are observed and cool core sys- 
tems are seen to have low temperatures. Defined in this way, 
At is positive for systems with core temperatures cooler 
than isothermal and negative for systems with core temper- 
atures warmer than isothermal. 

The entropy cut we have chosen (S40 = 50 keVcm 2 ) for 
dividing cool cores from warm cores roughly translates to a 
cut in central temperature gradients of At= 0.05 for relaxed 
systems. 



3.2.3 apparent dynamical state 

Each projection for each system is also classified into in- 
stances where it would exhibit no substructure and regular 
isophots to an observer, significant substructure, or so much 
substructure as to render the system indescribable by a f3- 
model (generally due to the close proximity of the secondary 
or various other brief but extreme distortions) . We refer to 
these situations as undisturbed, disturbed and strongly dis- 
turbed respectively. Classification is performed through vi- 
sual inspection of simulated 50ks Chandra observations for 
each output (with the system at z = 0.1), in each projec- 
tion. Our simulated observations were produced following 
the procedure detailed in PooleOO. 

Since instances of strongly disturbed morphology can 
not be assigned meaningful core radii, they are omitted from 
this morphological discussion. Such instances are rare and 
their exclusion does not represent a significant omission, al- 
though it is important to keep in mind that this exceptional 
class does manifest occasionally at early stages of the inter- 
action. 



3.3 Three Dominant Morphologies 

In Fig. [T]we present a At — r c plot for our nine simulations. 
Projections in y and z are displayed for each of the 120 
outputs (spanning 12 Gyr in 100 Myr intervals) of our sim- 
ulations with instances of undisturbed and disturbed mor- 
phology indicated by filled circles or crosses respectively. Our 
cuts in r c and S40 separate four potential populations which 
we identify in this plot using colours: compact cool cores 
(CCCs) are shown in blue, compact warm cores (CWCs) 
are in red, extended cool cores are in magenta and extended 
w arm cores (EWC s ) are green. We also plot the catalogue 
of lSanderson et al.l (120061 ) in black for comparison. 

Fig. [l] shows that very few instances of extended cool 
cores occur. Looking carefully at those circumstances which 
have been identified as an extended cool core, we find that 
in many cases it is a result of one of two effects. In some 
instances when the system appears disturbed, the passage 
of a merging core near to the primary core (either during 
first or second pericentric passage) can enlarge estimates of 
r c while the primary core remains sufficiently undisturbed 
to be cool. The system is very obviously disturbed at these 
times but not so much that a /3- model can not be reasonably 
fit to its surface brightness. The only case in which signifi- 
cant numbers of relaxed extended cool core states occur is 
the 10:1 head-on merger. The system has relaxed to a state 
with core radii only slightly larger and central entropies only 
slightly less than our morphology cuts in this case. 

With exceptional instances such as these aside, we see 
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Figure 1. Central temperature gradient (At) plotted against core radius (r c ) for each output of our nine simulations in y and z- 
projections. The cuts in core radius we use to separate compact-core systems (r c < 100 kpc) from extended-core systems (r c > 100 kpc) 
as well as the cut in Ay which roughly separates cool core (S40 < 50keVcm 2 ) from warm core (S40 > 50keVcm 2 ) systems are shown 
with dashed lines. Compact core systems with low/high central entropies (i.e. CCC and CWC systems) are labeled in blue and red 
respectively while extended core systems with low/high central entropies (i.e. ECC and EWC systems) are in magenta and green 
respectively. The quadrants of this plane corresponding to each morphology are labelled in the central pannel. Each population is 
further subdivided into systems showing obvious evide nce of substructu r e (cro sses) and systems with regular isophots (solid circles). 
Data from the statistically representative catalogue of ISanderson et alj 1120061 ) is shown in black with error bars (closed circles and 
crosses indicate systems noted in their study to be relaxed or visibly disturbed respectively). Black text around the boundary indicate 
the mass ratio and vt/V c of the simulations presented in each panel. 
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Figure 2. Central cooling time (t coo i) plotted against central temperature decrement (Ay) for each output of our nine simulations 
in y and z-projections (following the same format as Fig. [TJ. The cuts in Ay and t coo i which roughly separate cool core systems 
from warm core systems are shown with dashed li nes (the quadrants corr esponding to warm and cool morphologies are labelled in the 
central panel) while the i coo ; = 5 Gyr cut used bv lSanderson et al.l 1120061) to divide short cooling time systems from long cooling time 
systems is shown with dotted lines. Compact core systems with low (S40 < 50keVcm 2 ) and high (S40 > 50keVcm 2 ) central entropies 
(i.e. CCC and CWC systems) are labeled in blue and red respectively while extended core systems with low/high central entropies 
(i.e. ECC and EWC systems) are in magenta and green respectively. Each population is further subdivided into systems showing 
obvious evidence of substructure (crosses) and systems with regular isophots (solid circles). Data from the statistically representative 
catalogue of Sanderson et al. (2006) is shown in black with error bars (closed circles and crosses indicate systems noted in their study 
to be relaxed or visibly disturbed respectively). Black text around the boundary indicate the mass ratio and vt/V c of the simulations 
presented in each panel. 
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that throughout the duration of a merger an observer would 
classify the system as existing predominantly in one of three 
states: compact systems with cool or warm cores (CCC or 
CWC) or extended warm cores (EWC). In each case, the 
system may or may not also possess obvious substructure 
and/or disturbed isophots, meaning that members of each 
class could be found which appear to be disturbed or undis- 
turbed. Including the approximate cut in At which results 
from our cut in S40, we see that these three populations 
(especially the undisturbed set) map well to three separate 
regions in the Ar-r c plane. Encouragingly, we also find 
good agreement between the regions of this plane occupied 
by our simulations and by the observations. In Section [3.41 
we will examine the evolution of our systems through this 
plane and in Section [3.51 examine the processes principally 
responsible for setting a system's morphology. 

It has been noted that clusters with warm cores can 
have cooling times less than their median formation age 
i|Donahue et al.ll2005l ; [Sanderson etalll2006h . In Fig. we 
plot the cen tral cooling time (meas ured at 40 kpc) against 
At from the lSanderson et al.l (120061 ) catalogue and for each 
of our simulations. We can see from this figure that our 
simulations capture the relationship betw een cooling time, 
core t emperature and core size observed bv lSanderson et al.l 
(2006): warm cores can have short cooling times if they are 
compact. This is a straightforward consequence of the strong 
dependence of cooling on density and its weak dependence 
on temperature. 

Although compact cores are not necessarily cool they 
nevertheless have short cooling times (even when disturbed) , 
preserving the importance of radiative cooling in establish- 
ing their structure and evolution during and after a merger. 
Furthermore, although extended core morphologies gener- 
ally (but not necessarily) have longer cooling times than 
compact core systems, only one of the cases we have studied 
produces states with cooling times greater than 5 Gyrs for 
any significant length of time: our 3:1 Vt/V c = 0.15 case. 



3.4 Temporal Evolution of Morphology 

To see how the morphologies of our simulations progress 
through these three states, we present the temporal evo- 
lution of r c and At in Figs. [3] and [4] In both cases, a 
black curve traces the system's evolution as seen in the z- 
projection with coloured points plotted (following the for- 
mat of Fig.[T} to indicate each interval's morphological state. 
These curves change very little when viewed in y-projections 
and differ only at early times in the x-projection due to the 
apparently close proximity of the two cluster cores prior to 
interaction (the system can appear as a ECC system in x- 
projections prior to t c i osest ). To clarify these plots, we show 
only the ^-projection. We also note tciosest and t re i ax for 
each simulation to indicate the interval over which the sys- 
tem will tend to look disturbed. 



pericentric passage. For the head-on mergers, we see that 
in the 1:1 and 3:1 cases this lasts (with varying amplitude) 
until roughly t accrete while in the 10:1 case it lasts for only 
~0.5 Gyr. These elevated values of r c are produced by the 
explosive expansion of the core after the secondary's impact. 

In off-axis cases, these elevated values are primarily a 
product of the close proximity of the cluster cores during 
pericentric passage. In these cases there is an interval be- 
tween tdoscst and t a ccrete during which the system returns 
to a compact core state. This indicates that the initial in- 
teraction of the cores fails to disrupt the primary system to 
an extended core state in any of our off-axis mergers. 

Following taccrete when the remainder of the secondary 
core accretes to the remnant core, many of our merger rem- 
nants are in, or nearly in, extended core states. They sub- 
sequently evolve back towards their initial core size (rapidly 
in some cases) with only the 3:1 head-on and 1:1 and 3:1 
vt/Vc — 0.15 cases spending prolonged periods of time in 
extended core states. Of these, only the 3:1 v t /V c = 0.15 
case fails to formally recover a compact core by the end of 
the simulation. The remnant in this case reaches an extended 
core state with r c ~ 140 kpc following taccrete and evolves 
very little afterward. 

In Section 13.51 we will discuss the physical processes 
which shape the diversity in remnant core states. 



3.4.1 



core size 



In Fig.[3]we consider the evolution of r c for our simulations. 
We see that in all cases the system starts with a compact 
core (by construction) and experiences a period of substan- 
tially increased core size (often r c > 300 kpc) during first 



3.4-2 central temperature gradients 

In FigUwe plot the evolution of At for our simulations fol- 
lowing the format of Fig. [3] In all cases our systems begin 
with strong positive central temperature gradients (by con- 
struction) which remains intact until t , modulo a very slow 
decline due to the numerical effects discussed in Section [274] 

Considering the head-on mergers first, we see that in 
the rare 1:1 head-on case, At remains above at nearly 
all times. In the more common 3:1 and 10:1 head-on cases, 
the system briefly develops a negative central temperature 
gradient after tciosest and oscillates with a warm core two or 
three times on dynamical timescales. 

In off-axis cases, the primary core's positive central tem- 
perature gradient survives until taccrete- This indicates that 
the first interaction of the cores in off-axis cases does not 
heat the primary core to a warm state. At (or shortly fol- 
lowing) taccrete however, all of our 1:1 and 3:1 mergers expe- 
rience a rapid decline in At- In the 10:1 cases there is little 
evolution in At following taccrete due to the substantial dis- 
ruption of the secondary system prior to second pericentric 
passage in these cases. 

Hence, we find that cool cores resiliently maintain their 
integrity until a merging system's core makes direct impact, 
either at t c i OS est in head-on cases or during second pericen- 
tric passage at taccrete in 1:1 and 3:1 off-axis cases. Extended 
cores are not generically produced (except as very short lived 
transient states during which the system would appear visi- 
bly disturbed) and the cores generally recover cool states 
once relaxed (the most significant exceptions are the 3:1 
head-on and Vt/V c = 0.15 cases). 

We explore the physical processes which shape the di- 
versity of remnant core states in Section 13.51 
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Figure 3. The temporal evolution of core radius (r c ) for our simulations. Black curves illustrate r c (t) and coloured points indicate the 
morphology of the system measured in the ^-projection. Compact core systems with low (S40 < 50keVcm 2 ) and high (S40 > 50keVcm 2 ) 
central entropies (i.e. CCC and CWC systems) are labeled in blue and red respectively while extended core systems with low/high 
central entropies (i.e. ECC and EWC systems) are in magenta and green respectively. The cut in r c used to separate compact core 
systems from extended core systems is shown with dashed lines (see labels in central panel). Vertical dotted lines indicate (from left to 
right) t c i osest , taccrete and t re i ax for each case. The interval between t c i osest and t re i ax indicates the time during which the system 
will show signs of being disturbed. Black text around the boundary indicate the mass ratio and vt/V c of the simulations presented in 
each panel. 
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Figure 4. The temporal evolution of the central temperature gradient (Ay) following the same format as Fig. [3] Black curves 
illustrate Ar(t) and coloured points indicate the morphology of the system measured in the ^-projection. Compact core systems with 
low (540 < 50keVcm 2 ) and high (S40 > 50keVcm 2 ) central entropies (i.e. CCC and CWC systems) are labeled in blue and red 
respectively while extended core systems with low/high central entropies (i.e. ECC and EWC systems) are in magenta and green 
respectively. The Aj-= 0.05 line roughly separating cool core systems from warm core systems (classified by entropy) is shown with 
dashed lines (see labels in central panel). Vertical dotted lines indicate (from left to right) t c i oaest , taccrete and t re i ax for each case. The 
interval between t c i oaeat and t re i ax indicates the time during which the system will show signs of being disturbed. Black text around 
the boundary indicate the mass ratio and vt/Vc of the simulations presented in each panel. 
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3.5 The recovery of compact cool cores 

As we can see from both Figs. [3] and [4] a pattern emerges for 
the durations of EWC states (illustrated in green) produced 
in our mergers: near-axis (head-on and v t /V c — 0.15) and 
intermediate mass ratios (3:1) are the most efficient cases 
for producing this morphology. To understand this, we now 
seek to explain two things: how the initial conditions of the 
merger (namely, the mass ratio and orbits of the clusters) 
dictate the final remnant state which emerge at t re i ax and 
what properties of the remnant core's structure principally 
govern the recovery of CCC morphologies. 

In essence, this reduces to an examination of the final 
balance between shock heating and radiative cooling for ma- 
terial destined to end-up in the core of the remnant. Thus, 
the pattern of sustained EWC morphologies seen in Figs. [3] 
and 2] can be understood as the result of a "Goldie Locks" 
scenario where the circumstances necessary for effective core 
heating are optimized while the formation of dense struc- 
tures capable of offsetting heating through radiative cooling 
are prevented. 

These competing processes are both sensitive to the 
evolving dynamical state of the initial systems' cores and 
the relative densities and mach numbers of gaseous material 
forced to interact during the event. We find that a great va- 
riety of processes are involved in establishing this balance, 
with the effects of both shock heating and radiative cooling 
being a sensitive function of the mass ratio and initial orbits 
of the interacting clusters. 

For the remainder of this section we shall examine this 
variety. However, as discussed in Section [2U1 our sparse sam- 
pling of M p : M s and Vt/V c prevents us from systematically 
studying many of these processes in a quantitative way. We 
will thus restrict ourselves to a heuristic discussion in what 
follows. 



3. 5. 1 entropy and the formation of extended cores 

The temperature structure of a relaxed cluster's ICM is pri- 
marily set by its gravitational potential. Because of this and 
the strong p' 2 g dependence of the ICM's X-ray emissivity, 
variations in the central X-ray surface brightnesses of clus- 
ters (and hence, variations in r c ) are primarily a product 
of variations in central gas densities. Because it is directly 
tied to the mass accretion history of the system and the rel- 
ative balance between heating (from accretion shocks and 
feedback) and radiative cooling, these density variations are 
most straightforwardly interpreted in terms of variations in 
central (or minimum) gas entropies. 

In Fig. [5] we present the temporal evolution of the cen- 
tral entropies of our simulations where we see a good cor- 
relation between Saq{€) and r c (t) following treiax- Hence, 
the surface brightness core radii of our merger remnants are 
tightly related to their central entropy. We now address the 
question: what processes set this entropy? 

As in the plots of r c (t) and At(£), we observe important 
differences in 5*4o(t) between head-on and off-axis cases. In 
the head-on cases, one very dramatic increase in Sao occurs 
at t c i oses t followed by several small and erratic oscillations 
which end ~ 1 Gyr before t Te i ax . The details of this evolu- 
tion depend quite sensitively on the mass ratio of the event. 
Examining the dynamical evolution of these mergers in the 



spatial distribution of their surface brightness, temperature 
and entropy (Figs. 1 to 9 in Poole06) we can see why this is 
so. In the 3:1 and 10:1 cases, the secondary core fully pen- 
etrates the primary core triggering explosive shocks which 
heat the entire primary core from the inside out. In the 1:1 
case, the merging cores inelastically collide and remain in 
very dense structures following tdosest- This has the dual ef- 
fect of reducing the efficiency of entropy production through 
shocks and increasing the rate of radiative cooling. Together, 
these effects lead to a counter intuitive result: the 3:1 and 
10:1 head-on mergers ultimately heat the primary core more 
than the 1:1 head-on merger. 

In off-axis cases, there are slight and brief increases in 
central entropy at t c i osest in some cases, but there is gener- 
ally very little sustained increase in Sao until t accrete- This is 
because the initial interaction of the systems has very little 
effect on the primary core and is consistent with what we 
noted in Section[3]4] the primary core retains its initial CCC 
morphology until second pericentric passage in off-axis cases. 
In the 1:1 cases, the secondary core remains significantly in- 
tact until taccrete and we thus see a significant jump in en- 
tropy. In other words: the evolution of the 1:1 v t /V c = 0.15 
case following taccrete is very similar to the 3:1 head-on case 
following tciosest- In the 10:1 off-axis cases, the secondary 
core is so significantly disrupted by taccrete that essentially 
no evolution in central entropy is seen throughout the entire 
simulation. 

In the 3:1 off-axis cases, an evolution distinctly different 
from the 1:1 and 10:1 off-axis cases is observed. In Poole06 
we used our 3:1 Vt/V c = 0.15 case to illustrate how an off- 
axis merger of a cold-core system can generate low-entropy 
streams of material which accrete to the remnant core with 
high velocities. Fig. [S] illustrates that the accretion of these 
streams can have a significant impact on the central entropy 
of off-axis merger remnants. In both off-axis 3:1 cases we see 
a very low amplitude jump at taccrete when the surviving 
portion of the secondary core returns for its second pericen- 
tric passage. It is not until significantly after taccrete that 
peak levels of Sao are reached. These peaks occur roughly 
when the streams formed from the secondary cores are ac- 
creting most significantly onto the remnant. 

The heating effects of an accreting stream depend sen- 
sitively on its total mass, density, and mach number as well 
as the state of the remnant core onto which it falls. The ram 
pressure and tidal forces experienced by the cores of inter- 
acting systems, which are responsible for stripping the sec- 
ondary, forming the stream and disturbing the primary core, 
depend very sensitively on the initial mass ratio and impact 
parameter of the merger. The stream's mach number during 
accretion also depends on the amount of dynamical friction 
experienced by the secondary system during its orbit, and 
thus, the mass ratio of the interacting systems. Hence, the 
heating effects of stream accretion are very sensitive to the 
initial mass ratio and orbits of the interacting clusters. 

Thus, the central entropy of the remnant, which is a 
very sensitive function of the mass ratio and initial orbit 
of the merger, plays a key role in setting its core size. In 
the cases we have studied, the only long lived extended core 
remnant results from a slightly off-axis 3:1 merger for which 
stream heating is the dominant mode of core entropy pro- 
duction. 



12 Gregory B. Poole et al. 




2 4 6 B 10 



Figure 5. The temporal evolution of the central entropy (measured at 40kpc) following the same format as Fig. \3\ Black curves 
illustrate S4o(i) and coloured points indicate the morphology of the system measured in the z-projection. Compact core systems with 
low (S40 < 50keVcm 2 ) and high (540 > 50keVcm 2 ) central entropies (i.e. CCC and CWC systems) are labeled in blue and red 
respectively while extended core systems with low/high central entropies (i.e. ECC and EWC systems) are in magenta and green 
respectively. The cut in S40 used in our morphological classification to separate cool core systems from warm core systems is shown 
with dashed lines (see labels in central panel). Vertical dotted lines indicate (from left to right) t c i oseat , taccrete and t re [ ax for each 
case. The interval between t c i osest and t re i ax indicates the time during which the system will show signs of being disturbed. Black text 
around the boundary indicate the mass ratio and vt/V c of the simulations presented in each panel. 



3.5.2 cooling timescales and cool core recovery 

Radiative emission can reduce the temperature of material, 
but must it do so? Imagine a disturbed parcel of low-entropy 
gas radiatively cooling as it falls to the center of the system. 
As it cools, its entropy will fall on its cooling timescale, 
tcooi on ST. However, because of its dynamically dis- 
turbed state, its temperature can drop only if the cooling 
timescale is sufficiently short with respect to its local dynam- 
ical time, tdyn °c (GAp c ) -0,5 . If t coo i 3> tdyn, temperature 
increases from compression dominate over radiative losses. 
Under these circumstances, the temperature can only fall 
once the material reaches a quasi-static state. It will then 
do so on timescales governed by the cooling time. 

At all times following taccrete and for all the cases we 
have studied, we find that t coo i 3> tdyn within r c (by ap- 
proximately an order of magnitude) for our disturbed but 
relaxing remnants. As a result, cooling does not lower our 
remnant cores' temperatures until they have returned to 
quasi-stable states of hydrostatic equilibrium and a substan- 
tial fraction of a cooling time has passed. In Poole06 we 
examined how a system returns to hydrostatic equilibrium 
following a merger. We showed that our remnants gener- 
ally return to within 10% of hydrostatic equilibrium at -R500 
0.5 — 1.0 Gyrs following the system's apparent relaxation at 
treiax- Examining this behavior at -R2500 (roughly r c ) we find 



that the system returns to within 10% of hydrostatic equi- 
librium a little sooner and generally within 0.5 Gyrs of when 
it looks relaxed at treiax- Hence, cooling generally will not 
have a significant effect on central remnant temperatures 
until a time t CO oi (i.e. ~ 2 — 5 Gyrs) after appearing re- 
laxed at t re i ax . This accounts for the occurrence of compact 
warm core (CWC) morphologies as well as the timing of the 
late increases in At occurring in several of our 1:1 and 3:1 
merger remnants (whose central cooling times are short). 

To conclude, we find that following a disturbance from 
a merger, the initially dense cores of CCC systems gener- 
ally recover much faster (roughly a few dynamical times fol- 
lowing to) than their cool temperatures (roughly t coo i after- 
ward) . The existence of compact war m core (CWC) systems 
observed by Cha ndra and XMM (jSanderson et al.l 120061 ; 
IPratt et al.l 120071 ) can thus be understood naturally within 
the framework of hierarchical structure formation without 
appealing to other heating sources such as AGN. Their exis- 
tence is merely a natural consequence of hierarchical struc- 
ture formation and the fact that t CO oi 3> tdyn in cluster cores. 



3.6 Morphological segregation in scaling relations 

Scaling relations between mass, X-ray temperature and lu- 
minosity are powerful tools for studying the processes that 
shaped the structure of the ICM. It has been known for some 
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time (|Fabian et al.lll994J ; iMarkevitchll 19981 ) that a significant 
fraction of the scatter in these relations is due to variations 
in the core propert ies of clusters. Thi s has been supported by 
recent theoretical (JQ'Hara et alj2006L Poole07, McCarthy et 
al. 2007) and ob servational studies (jSanderson et al.l 120061 ; 
IPratt et a.1.1 120071 ). For this reason, and because the gener- 
ation of accurate mass functions from scaling relations has 
been a significant motivator for their study, past analyzes 
have typically been conducted with the central cores of sys- 
tems excised to reduce scatter. 

However, it has been shown by M04 that there is inter- 
esting structure in these relations which correlate with the 
core properties of clusters. They use the ACC catalogue of 
iHorner I l|200ll ) to compare their analytic entropy injection 
model to published Chandra and XMM temperature pro- 
files, classifying them by the presence or absence of central 
temperature gradients. They further separate these systems 
into classes showing evidence of significant or little substruc- 
ture. Despite the small number of published Chandra and 
XMM temperature profiles available to them at the time, 
they found that these classes clearly (but not necessarily 
distinctly) segregate in both the L x — T x and /3 — r c planes. 

In Poole07 we examined the evolving global properties 
(including X-ray temperature and luminosity) of our simu- 
lations, as well as the scaling relations generated from them. 
Having quantified the evolving morphologies of our cluster 
cores, we now examine how the three morphology classes 
we have identified segregate in these planes. In Poole07 we 
found that the evolution in the L x —T x , Mt — L x , Mt — T x , 
SZ — L x and SZ — T x planes exhibited a generic pattern 
and we find the same for morphological segregation in these 
planes as well. For this reason we will only examine the 
L x — T x and (3 — r c relations here. 



3. 6. 1 the L x — T x plane 

In Fig. [6] we present the positions of our systems in the 
L x —T x plane, at each time in y and z-projections, when they 
appear as CCC (blue), CWC (red) or EWC (green) systems. 
The initial and final states of the systems are indicated with 
black four and five point stars respective ly. These are co m- 
pared against the observed catalogue of IHorner I (|200ll . in 
black) fo r which we use the cata l og of R OS AT 1-D /3-model 
fits from iReiprich fc Bohringer I ((2002) to identify systems 
with compact cores (r c < 100 kpc) using filled circles and 
systems with extended cores (r c > 100 kpc) using open cir- 
cles. 

This plot reinforces what was found in Poole07: the low- 
luminosity side of the L x — T x scaling relation can not be 
significantly populated with the mergers we have studied. 
It is comforting to note that our systems remain confined 
within the scatter of the observations at all times, but in- 
stances of systems lying near the low luminosity side of this 
scatter are rare and short lived. 

Comparing the regions occupied by systems which ex- 
hibit compact and extended cores, we find a good agreement 
between the observations and our simulations (with the ex- 
ception of the rare 1:1 off-axis mergers, which yield much 
more under- luminous compact core systems). Furthermore, 
there is a very clean divide between warm and cool cores 
in this plane. We have plotted a fiducial entropy injection 
model with a minimum entropy of S — 150keVcm 2 (com- 



puted following the general approach of McCarthy et al. I 
12004 . see Poole07 for some additional details) to roughly 
quantify the location of this divide. 

For our simulations, instances of disturbed morphology 
generally lie near or above this fiducial line while relaxed 
morphologies lie below it. This is consistent with common 
anecdotal observations that extended core systems gener- 
ally appear disturbed. However, our results suggest that dis- 
turbed morphologies may correlate better with the presence 
of warm cores than with the presence of extended cores. It 
is important to emphasize however that our simulations do 
not populate the low- luminosity side of the scatter in this 
plane nor do they yield relaxed extended core systems such 
as Abell 2034 or Abell 2631. Trends such as these suggest 
that the observed scatter in the L x — T x plane is a result of 
physical processes other than mergers. 

We have similarly examined the other scaling relations 
discussed in Poole07 and find the same qualitative behavior 
in all cases. 



3.6.2 the/3- 



plane 



In Fig. [7] we present the positions our simulations occupy in 
the P — r c plane, at each time in y and z-projections, when 
they appear as CCC (blue), CWC (red), ECC (magenta) or 
EWC (green) systems. Instances of relaxed-looking morphol- 
ogy are marked with solid circles and disturbed morphol- 
ogy with crosses. The initial and final states of the systems 
are indicated with black four and five point stars respec- 
ti vely. These are compared agai nst the observed catalogue 
of IReiprich fc Bohringerl (|2002l ) (black dots). We identify 
systems with short cooling times (t coo i < 5 Gyrs) with large 
filled circles and systems with long cooling times (t coo i > 
5 Gyrs) with large open circles. These cooling times a re ob- 
tained by combining t emperatures from Horner I (120011 ) with 
denisties derived from IReiprich fc Bohringer I ( 2002j) . 

We can see from this plot that /3 and r c express a clear 
correlation in the observations and that our mergers repro- 
duce some of its features. In particular, our simulations ac- 
count for the high-/3 envelope of the observed scatter in this 
plane. Furthermore, our simulations populate the narrow 
region occupied by the observations at the largest values of 
r c . While they significantly populate this region, they are al- 
ways visibly disturbed. Observed systems in this region for 
which we have dynamical classifications (|McCarthv et al. I 
|2004J) are similarly all disturbed. 

However, our simulations fail to produce states which 
last for significant lengths of time in the region of this plane 
occupied by the vast majority of observed extended core 
systems. Most observed extended core systems have values of 
P substantially less than the states produced by our mergers. 

In Fig. [2] we show that our mergers rarely produce sys- 
tems with cooling times in excess of 5 Gyrs. This is primarily 
a consequence of the cores of our systems remaining compact 
for most of their evolution. We can see in this plane that ob- 
served systems with short cooling times are generally found 
along the high-/? side of the observed scatter in this plane, 
in good accord with our simulations. Observed systems with 
tcooi > 5 Gyrs lie in the low-/3/large-r c region of the plane 
which our simulations fail to significantly populate. 

To conclude, our simulations are able to account for the 
high-/? envelope occupied by observed systems with short 
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Figure 6. A plot of L x against T x illustrating the morphological segregation in this plane. Compact core systems with low (S40 < 
50keVcm 2 ) and high (S40 > 50keVcm 2 ) central entropies (i.e. CCC and CWC systems) are labeled in blue and red respectively while 
extended core systems with high central entropies (i.e. EWC systems) are in green. Filled circles and crosses mark syste ms which 
appea r relaxed or unrelaxed (respectively) in simulated 50 ks Chandra observations. Black dots are the observed catalogue of iHorner I 
I l200ll) with open circles marking systems with extended cores (r c > 100 kpc) and filled circles marking systems with compact cores 
(r c < 100 kpc). The black dashed line is a fiducial entropy injection model with So = 150keVcm 2 (see Poole07 for details). Black text 
around the boundary indicate the mass ratio and vt/V c depicted by each panel. 
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cooling times as well as the tail of disturbed systems at the 
highest values of r c . They do not account for the bulk of 
observed extended core systems with long cooling times. 



4 METALLICITY GRADIENTS 

It is a well established fact that extended core systems 
tend to lack metallicity gradients (with mean amplitudes 
of 0.1 — 0.2^0 ) while CCC systems routinely exhibit cen- 
trally peaked metallicity gradients within 0.27?igo, often 
reaching central values whi ch approach solar metallicity 
|De Grandi fc MolendU 120021 ). This, along with anecdotal 
evidence suggesting that extended core systems tend to be 
preferentially disturbed has been used to claim, as a result, 
that mergers effectively mix the ICM. 

To examine this hypothesis, we have initialized our 
systems to have the meta llicity profiles presented by 
|Pe Grandi fc Molendil (J2002J) and examined the resulting 
metallicity profiles of our merger remnants. We present these 
results in Fig. [5] where the initial metallicity profiles of our 
primary systems are presented in blue and the final metal- 
licity profiles of our remnants in red. 

It is clear from this figure that mergers do not efficiently 
mix the material of the initial systems. This makes sense 
given the fact that the enhanced metallicity is initially tied 
to the lowest entropy material in the clusters. This material 
is not preferentially heated above the entropy of the sur- 
rounding material and thus, remains the low entropy mate- 
rial constituting the remnant. Once the system relaxes, the 
low-entropy/high-metallicity material sinks to the center of 
the system due to a lack of buoyant support, reestablishing a 
central metallicity peak with only a slightly lower amplitude 
and slightly broader distribution. 



5 DISCUSSION 

In this paper, we have shown that the notions of compact- 
ness, coolness, and active cooling for cluster cores are related 
but not as si mply as once thou g ht. A s illustrated by the ob- 
servations of ISanderson et al.l i|2006h , compact cores need 
not be cool and can have central cooling times as long as 
4Gyrs while warm cores can have cooling times as short as 
2 Gyrs. We have shown how several aspects of these trends 
can be accounted for naturally in theories of hierarchical 
structure formation: radiative cooling naturally produces re- 
laxed clusters with compact cool cores and significant merg- 
ers generally disturb them into warm core states with a 
range of core radii. The central entropies generated by a 
merger event principally determine the core size of the rem- 
nant and if the resulting cooling times are short, the system 
will evolve back towards a CCC state after appearing re- 
laxed. The core dynamically relaxes much more quickly than 
it radiatively cools so a compact core recovers its central den- 
sities (and hence, surface brightness core size) quicker than 
it's cool temperatures return. This naturally explains the 
population of compact warm cores observed by Chandra. 

Our mergers do not produce relaxed clusters with large- 
r/low-/? X-ray surface brightness profiles and long cool- 
ing times. They do not erase metallicity gradients either. 
These results imply that once established, compact cool 



cores are remarkably resilient to disturbances from mergers. 
How then are relaxed systems with high central entropies 
(S40 > 150keVcm 2 ) and extended surface brightness cores 
(r c > 150 kpc) produced? In Poole07 we suggested the pos- 
sibility that mergers could induce AGN activity, augmenting 
the heating from a m erger event. Significant o utbursts such 
as those observed by iMcNamara et al.l (|2005r i could poten- 
tially heat the cores of a disturbed system to EWC states 
but even MS0735.6+7421, the largest AGN outburst yet ob- 
served, still has a significant metallicity gradient. It would 
appear that neither AGN outbursts nor mergers are able to 
sufficiently erase a metallicity gradient once it is established. 

How then do we reconcile the lack of metallicity gradi- 
ents in observed EWC systems with our results? It may be 
that mergers yield more mixing than our simulations sug- 
gest. SPH notoriously handles turbulence poorly which may 
be leading us to underestimate the degree of mixing in merg- 
ers. On the other hand, the non-diffusive nature of SPH pre- 
vents sedimentation which would otherwise reenforce metal- 
licity gradients. 

Alternatively, we may need to seek mechanisms by 
which cool cores are prevented from forming in the first 
place. One means may be through variations in mass ac- 
cretion histories. We find that major near-axis mergers heat 
cores most efficiently and at early epochs, such events (al- 
though between much lower mass systems tha n we study) 
are ex pected to be more common. Recently iBurns et al.l 
(|2007l ) have argued that early major mergers prevent com- 
pact cool cores from forming. Systems having experienced a 
greater proportion of such events may have altered star for- 
mation histories as a result, yielding the observed differences 
in metallicity profiles between compact a nd extended core 
system s. On the other hand, as noted by iMcCarthv et al.l 
(2007), the virtual absence of extended core systems in 
all other studies involving cosmological simulations is com- 
pelling support for the idea that mergers alone are not 
responsible for the full range of observed cluster states. 
Clearly, further investigation is required to enable us to 
select between competing possibilities. Towards this goal, 
we should keep a keen eye-out for relaxed EWC systems 
with t coo i > 5 Gyrs. Such a system would be a strong can- 
didate for a cluster that has been significantly heated by 
non-gravitational processes and would shed a great deal of 
light on these issues. 



6 SUMMARY AND CONCLUSIONS 



We use the simulations presented in IPoole et al.l 1)20061) to 
examine the effects of mergers on the properties of cool cores 
in X-ray clusters. Our results suggest a natural scheme for 
classifying the morphology of clusters based on their cen- 
tral surface brightness and entropy profiles and we show 
that this scheme effectively captures the diversity of clus- 
ter characteristics revealed by recent high-resolution X-ray 
observations. Our study illustrates how these morphologies 
naturally emerge in the context of hierarchical clustering. 
For our present purposes we have not included the effects of 
pre-heating or post-heating by AGN. We find: 

• Systems exist in three dominant morphologies: those 
with compact cores and central temperatures which are 
cool (CCC systems) or warm (CWC systems) and extended 
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Figure 8. A plot comparing the initial and final metallicity gradients of our simulations. Initial profiles of the primary are shown 
with thick red lines and the final remnant profile is shown with a thick black line. The final distribution of material originating in 
the secondary is sh own with a thi n blue line. In all cases, the metallicity profiles of our systems are initialized to have the metallicity 
profiles observed bv lDe Grandi &; Mo lcndi (20Q3). 



warm core (EWC) systems. Members of each population can 
appear either relaxed or disturbed. 

• Examining the temporal evolution of the morphology 
of our systems, we find that during merger events, CCC 
systems are only disturbed by direct collisions with incoming 
cores. This can occur in head on collisions or during second 
pericentric passage in off-axis mergers. 

• Our merger remnants are generally heated to warm core 
states by the time they relax. In cases where they relax to 
compact core structures, they typically recover cool cores 
(i.e. return to CCC states) after ~ 3 Gyrs. The central 
entropy of a merger remnant principally determines whether 
the system initially relaxes to a compact or extended core 
state. 

• A number of highly non-linear, strongly coupled pro- 
cesses are responsible for setting the central entropy of 
the remnant. The outcome is very sensitive to the initial 
mass ratio and orbits of the interacting systems. A heuristic 
overview of our simulation results are as follows: 

- In on-axis cases, 3:1 and 10:1 mergers paradoxically 
result in a more heated remnant core than the 1:1 mergers. 
This is due to enhanced heating from the penetration of 
the secondary core through the primary core and reduced 
cooling efficiencies in the low-density transient structures 
which result. In 1:1 head-on cases the cores inelastically 
collide and remain in high density structures. This re- 
duces the efficiency of entropy production by shocks and 
increases entropy loss through radiative cooling. 

- Of the off-axis mergers we have studied, the case 



which generates the most significantly heated remnant 
core is a near-axis 3:1 merger. In this case, the major- 
ity of core heating is due to the accretion of a low-entropy 
stream of material formed from the disruption of the sec- 
ondary core. 

• Since td yn "C t coo i for all our relaxing merger remnant 
cores, compression prevents their core temperatures from 
falling until after they relax to the compact states allowed 
by their remnant central entropies. This naturally explains 
the population of relaxed CWC systems observed in recent 
Chandra and XMM observations without any need to in- 
voke AGN feedback. We do not, however, exclude the pos- 
sibility that extreme AGN outbursts may on rare occasions 
transform compact cold cores into compact warm cores. 

• The morphological s e gregat ion in scaling relations 
noted by iMcCarthv et al. I (J200J) is qualitatively reflected 
in the results of our mergers: EWC systems preferentially 
lie towards the low-luminosity/high-temperature side of the 
scatter in the L x — T x plane. However, as noted in Poole07, 
the typical cluster mergers we have studied (3:1 and 10:1 
cases) do not produce systems with sufficiently high cen- 
tral entropies to account for the most under-luminous EWC 
systems. 

• In the p — r c plane, several trends in the observations 
can be naturally accounted for by mergers, but some can 
not. In particular, our mergers nicely account for the ob- 
served dependence on r c of the high-value limit of /3, with 
CCC systems lying at low values of r c and disturbed EWC 
systems lying at high values of r c . They also account for the 
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observed tail of disturbed systems at the largest observed 
values of j3. They can not account for extended core sys- 
tems with low-/3 and long cooling times (i.e. the majority 
of extended core systems) . 

• Mergers do not efficiently mix the ICM. Mergers be- 
tween systems which initially host central metallicity gra- 
dients produce merger remnants with metallicity gradients 
only slightly reduced in central amplitude and slope than 
the initial systems. 

Our results thus pose a significant challenge to a long 
held hypothesis in the study of X-ray clusters: that extended 
core systems are the remnants of mergers between relaxed 
compact core systems. They suggest that compact cool core 
systems are remarkably stable to disturbances from cluster 
mergers. Consequently, to explain the observed spectrum 
of X-ray cluster morphologies, further study of mechanisms 
capable of preventing the formation of compact cool cores is 
needed. 
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